
Calhoun 

iniQiuiic^iul Ar{hiv« of tilt Mil vdl Poii^roduiit School 


Calhoun: The NPS Institutional Archive 
□Space Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


2014-06 

Optimal sensor placement in active 
multistatic sonar networks 

Kuhn, Tobias Uwe 

Monterey, California: Naval Postgraduate School 


http://hdl.handle.net/10945/42665 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


htt p://w ww. n ps. e du/l ib ra ry 


Callwuo is the Naval Postgraduate School's public access digital repository for 
research mate rials and institutiional publicatkios created by the NPS community. 
Calhoun is named for Professor of Mathematics Guy K. Caftiouo, NPS's first 
appointed — and published — schoteily author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 Univefsity Circle 
Monterey, California USA 93943 







NAVAL 

POSTGRADUATE 

SCHOOL 


MONTEREY, CALIEORNIA 

THESIS 


OPTIMAL SENSOR PLACEMENT IN ACTIVE 
MULTISTATIC SONAR NETWORKS 

by 

Tobias Uwe Kuhn 
June 2014 

Thesis Co-Advisors: Emily M. Craparo 

Craig W. Rasmussen 

Second Reader: Miimtaz Karata§ 


Approved for public release; distribution is unlimited 




THIS PAGE INTENTIONALLY LEET BLANK 



REPORT DOCUMENTATION PAGE 

Form Approved 0MB No. 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instruction, 
searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments 
regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden to Washington 
headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202—4302, and 
to the Office of Management and Budget, Paperwork Reduction Project (0704-0188) Washington DC 20503. 

1. AGENCY USE ONLY (Leave Blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

06-20-2014 Master’s Thesis 11-06-2013 to 06-20-2014 

4. TITLE AND SUBTITLE 

OPTIMAL SENSOR PLACEMENT IN ACTIVE MULTISTATIC SONAR NETWORKS 

5. FUNDING NUMBERS 

6. AUTHOR(S) 

Tobias Uwe Kuhn 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Naval Postgraduate School 

Monterey, CA 93943 

8. PERFORMING ORGANIZATION REPORT 
NUMBER 

9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

N/A 

10. SPONSORING / MONITORING 

AGENCY REPORT NUMBER 

11. SUPPLEMENTARY NOTES 

The views expressed in this document are those of the author and do not reflect the official policy or position of the Department of 
Defense or the U.S. Government. IRB Protocol Number: N/A. 

12a. DISTRIBUTION / AVAILABILITY STATEMENT 

Approved for public release; distribution is unlimited 

12b. DISTRIBUTION CODE 

13. ABSTRACT (maximum 200 words) 

Recently the idea of deploying non-collocated sources and receivers in multistatic sensor networks (MSNs) has emerged as a promis¬ 
ing area of opportunity in sonar systems. This thesis addresses point coverage sensing problems in MSNs, where a number of points of 
interest have to be monitored in order to protect them from hostile underwater assets. We consider discrete “cookie cutter” sensors as 
well as various diffuse sensor models. By showing that the convex hull spanned by the targets is guaranteed to contain optimal sensor 
positions, we are able to limit the solution space. Using a cookie cutter sensor model, we are able to exclude even more suboptimal 
solutions by determining range-of-the-day, source and receiver circles. 

To address the nonconvex single-source placement problem, we develop the Divide Best Sector (DiBS) algorithm, which quickly 
provides an optimal source position assuming fixed receivers. Starting with a basic implementation of DiBS, we show how incorpo¬ 
rating advanced sector splitting methods and termination conditions further improve the algorithm. We also discuss two ways to use 
DiBS to find multiple source positions by placing sensors iteratively or simultaneously. Finally, we conclude that DiBS is a fast and 
simple algorithm that supports a wide variety of sensor models, various termination conditions, and objective functions. 

14. SUBJECT TERMS 

Multistatic Sensor Network, Point Coverage Sensing, Optimization, Divide Best Sector 

15. NUMBER OF 

PAGES 73 

16. PRICE CODE 

17. SECURITY 

CLASSIFICATION OF 

REPORT 

Unclassified 

18. SECURITY 

CLASSIFICATION OF THIS 

PAGE 

Unclassified 

19. SECURITY 

CLASSIFICATION OF 

ABSTRACT 

Unclassified 

20. LIMITATION OF 
ABSTRACT 

uu 


NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 


Prescribed by ANSI Std. 239-18 


1 




























THIS PAGE INTENTIONALLY LEET BLANK 


11 



Approved for public release; distribution is unlimited 


OPTIMAL SENSOR PLACEMENT IN ACTIVE MULTISTATIC SONAR 

NETWORKS 

Tobias Uwe Kuhn 
Major, German Army 

M.S., University of the German Armed Forces Munich, 2005 

Submitted in partial fulfillment of the 
requirements for the degree of 

MASTER OF SCIENCE IN OPERATIONS RESEARCH 
AND MASTER OF SCIENCE IN APPLIED MATHEMATICS 

from the 

NAVAL POSTGRADUATE SCHOOL 
June 2014 


Author: Tobias Uwe Kuhn 

Approved by: Emily M. Craparo 

Thesis Co-Advisor 

Craig W. Rasmussen 
Thesis Co-Advisor 

Miimtaz Karata§ 

Second Reader 

Robert Dell 

Chair, Department of Operations Research 
Carlos Borges 

Chair, Department of Applied Mathematics 



THIS PAGE INTENTIONALLY LEET BLANK 


IV 



ABSTRACT 


Recently the idea of deploying non-collocated sources and receivers in multistatic sensor 
networks (MSNs) has emerged as a promising area of opportunity in sonar systems. This 
thesis addresses point coverage sensing problems in MSNs, where a number of points of 
interest have to be monitored in order to protect them from hostile underwater assets. We 
consider discrete “cookie cutter” sensors as well as various diffuse sensor models. By 
showing that the convex hull spanned by the targets is guaranteed to contain optimal sensor 
positions, we are able to limit the solution space. Using a cookie cutter sensor model, 
we are able to exclude even more suboptimal solutions by determining range-of-the-day, 
source and receiver circles. 

To address the nonconvex single-source placement problem, we develop the Divide 
Best Sector (DiBS) algorithm, which quickly provides an optimal source position assuming 
fixed receivers. Starting with a basic implementation of DiBS, we show how incorporating 
advanced sector splitting methods and termination conditions further improve the algo¬ 
rithm. We also discuss two ways to use DiBS to find multiple source positions by placing 
sensors iteratively or simultaneously. Finally, we conclude that DiBS is a fast and simple 
algorithm that supports a wide variety of sensor models, various termination conditions, 
and objective functions. 
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Executive Summary 


Active sonar systems have long constituted an important sensing mechanism aboard sub¬ 
marines and ships in anti-submarine warfare (ASW). Recently, however, the idea of deploy¬ 
ing non-collocated sources and receivers has emerged as a promising area of opportunity 
in sonar systems. A multistatic sensor network (MSN) consisting of a number of non- 
collocated sources and receivers carries a number of advantages such as more complicated 
countermeasure tactics, improved deployment opportunities as well as reduced costs and 
noise pollution. The complex geometry of MSNs, however, makes deployment decisions 
like number of sensors and their locations more challenging. 

Out of three sensing problems studied in the literature (area search, barrier search, and 
point coverage) this thesis addresses point coverage problems assuming cookie cutter sen¬ 
sor models as well as diffuse sensor models. We further assume stationary targets and sen¬ 
sors in a two-dimensional Euclidean plane with homogeneous environment. Additionally 
we neglect the target’s aspect dependence and assume independence for all probabilities. 

We start by showing that for each sensor location outside the convex hull encasing all 
targets, there exists a position inside whose detection probability is at least as good. Based 
on this, we limit the search for optimal sensor locations to the convex hull of the targets. 

Next, we introduce the notion of range of the day circles (RDCs) for cookie cutter sensor 
models, i.e., circles around targets whose radii are the range of the day. Utilizing common 
points, we establish the minimal set of clusters of RDCs, G, whose cardinality represents 
a lower bound on the number of sensors required to cover all targets. Moreover, G lets us 
define an upper bound on the number of detected targets with a restricted number of sensors. 
Receiver circles (RCs) and source circles (SCs) provide additional means to narrow down 
possible sensor locations. 

Up to this point, we can formulate integer nonlinear programs (INLPs) with reasonable run 
time only to find optimal locations for exactly one source and one receiver, assuming a 
cookie cutter sensor model. Hence we develop the Divide Best Sector (DiBS) algorithm to 
find a single optimal source position assuming fixed target and receiver positions regardless 
of the applied sensor model. The algorithm divides the area of possible source locations 


XV 




into rectangular sectors and evaluates their upper bounds. At each iteration we divide the 
sector with the highest upper bound into smaller sector. Eventually, the algorithm reaches 
a termination condition such as maximum size of a sector’s longest edge or maximum 
optimality gap. 

We proceed to investigate and assess some details of DiBS starting with discovering that a 
small number of newly created sectors per iteration as well as fewer initial sectors greatly 
reduce the number of total sectors created and evaluated by the algorithm. This subse¬ 
quently leads to better run times. 

Another method to further improve the algorithm is rotating the plane such that the edges 
of the minimum-area rectangle encasing all targets are parallel to .r- and y- axis. By doing 
so, we minimize the sectors’ area outside the convex hull. 

Furthermore, we observe that uneven rectangles are able to increase the total number of 
sectors and therefore run time. We discuss two ways to address this issue by adjusting the 
termination condition and by creating only square sectors. 

Finally, we explore two approaches to apply DiBS to problems that require finding an 
optimal placement for multiple sources. Here, the iterative method provides a lower and 
upper bound assuming a cookie cutter sensor model. The second approach places sources 
simultaneously and evaluates all possible combinations of sectors. This approach, however, 
can only be applied to small numbers of sources. 
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CHAPTER 1: 
Introduction 


1.1 Background 

Active sonar systems have long constituted an important sensing mechanism aboard sub¬ 
marines and ships in anti-submarine warfare (ASW). In a typical sonar system, a ping is 
sent out and the echo yields information about other objects in the area. Recently, however, 
the idea of deploying non-collocated sources and receivers has emerged as a promising area 
of opportunity in sonar systems. 

A multistatic sensor network (MSN) consisting of a number of non-collocated sources and 
receivers carries a number of advantages. Cox (1989, p. 23) states that “countermeasure 
tactics are greatly complicated if the target does not know the position of the receivers.” 
This is justified by the fact that receivers do not send out pings and thus do not reveal 
their locations. Beyond that, procurement estimates of the United States Navy (USN) in¬ 
dicate that sources might cost about five times as much as receivers (USN, 2014). Thus, 
deploying more receivers than sources might significantly reduce costs without sacrific¬ 
ing performance. Also, a “multistatic system can employ different platforms for sources 
and receivers. A ship might be the source, while the receivers are sonobuoys” (Washburn, 
2010, p. 1). On top of this, a ping might be received by multiple receivers. Coon (1997) 
and Simakov (2008) discuss how to merge multiple detections into a single alert that is 
more precise and eliminates some of the false alarms that occur on traditional (monostatic) 
sonar systems. 

The performance of a bistatic sonar system, i.e., an MSN with exactly one source and one 
receiver, is significantly more difficult to model than a monostatic sonar system, leading 
to challenges in optimal deployment and usage. The primary source of these challenges 
lies in the differences of the geometry of both systems. In a monostatic sonar system the 
detection probability is mainly related to the distance between the sonar device and the po¬ 
tential target. This relationship is more complicated in a bistatic model, where the detection 
probability depends on the product of the target-source and target-receiver distances (Cox, 
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receiver r 


source i 



target t 


Figure 1.1: Geometry of a Bistatic Sonar System - The detection probability for target t 
depends on the product of the target-source and target-receiver distance. 

1989) as displayed in Figure 1.1. Moreover, the analytical challenges are exacerbated by 
MSNs involving multiple sources and receivers. 

1.2 Literature Review 

Researching sonar technology and MSNs in particular is a wide field with most diverse 
subcategories. We will discuss a selection in this section, starting with the geometry of 
multistatic models and following with examples for usage and deployment. 

1.2.1 The Geometry of Multistatic Sensor Networks 

Multiple authors describe the geometry of MSNs. Cox (1989) analyzes the relationship be¬ 
tween monostatic and bistatic active sonars. He derives that a detection probability contour, 
i.e., a contour consisting of all locations for a target t with the same detection probability, 
is defined by the constant product 



( 1 . 1 ) 


where pt,s,r is constant, and dt^s and dt^r are the distances from a target t to a source 5 
and to a receiver r respectively. The distances are also illustrated in Figure 1.1. Those 
contours are geometric figures known as Cassini ovals shown in Figure 1.2. The team of 
Matthew Fewell and Sylvia Ozols from the Australian Defence Science and Technology 
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^source; A receiver 



Figure 1.2: Examples of Cassini Ovals - Cassini ovals are shown with respect to source and 
receiver distance ds.r where po denotes the range of the day. If source and receiver are collocated 
then the Cassini oval is a circle that corresponds to a monostatic sensor. Increasing the distance 
changes the shape from oval over dog-bone shape up to two separate egg shapes. 


Organisation (DSTO) published multiple reports on MSNs. Fewell and Ozols (2011, p. 4) 
state that the constant pt,s,r from Equation (1.1) is the equivalent range to the range from a 
monostatic sonar system, which enables us to compute detection probabilities for bistatic 
sensor networks. The authors point out some issues with this model, e.g., “that not only [the 
detection probability is] high at a receiver but also the value is unaffected by the source- 
receiver distance. This seems inadequate (Fewell & Ozols, 2011, p. 7).” They list various 
possibilities to modify the model to account for these issues. 

There are multiple ways to determine Pt,s,r, the probability to detect target t with source s 
and receiver r in a multistatic model as displayed in Figure 1.3. The simplest one is the 
cookie cutter sensor model, also known as the definite range sensor in literature. Based on 
the range of the day po, i.e., the distance from a monostatic sensor to a target where the 
detection probability is 50%, we define Pt^s,r as 

_ Jl if pb < Pt,s,r < Po, 

[0 otherwise. 

Other names for the range of the day are the detection range or i? 5 o. The blind zone pt, 
is the area where targets cannot be detected because their echoes arrive at nearly the same 
time as the ping from the source. Then the total detection probability for target t, Pt, i.e.. 
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Figure 1.3: Sensor Models - The probability curves for the three sensor models cookie cutter, 
Fermi and exponential function are displayed. Range is expressed as multiples of the range of 
the day po. All models have a Pt^s,r = 0.5 at Pt,s,r/po = 1 to be consistent with the definition of 
the range of the day. 

the probability to detect target t with all sensors in scene, is defined as 

Pt= max Pt^s,r- (1-3) 

{s,r)eSxR 

This means, if at least one pair of sensors detects target t then /f = 1, otherwise Pt = 0. 
Although Equations (1.2) and (1.3) form an analytically convenient model, it lacks some 
features of real sensors such as a gradually decreasing Pt^s,r with increasing p. 

Hence, the DSTO team proposes two diffuse sensor models in (Fewell & Ozols, 2011). 
The first is the Fermi function^ 




{ _^_ 

1 ° 


if Pt,s,r ^ Pb^ 

Otherwise. 


(1.4) 


The diffusivity parameter b determines how rapidly probability values change when chang¬ 
ing range p. As Z? —)■ 0, the Fermi function approaches the cookie cutter model. 


^Discovered in 1926 by Enrico Fermi and Paul Dirac when researching electron behavior. 
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The second model is the exponential function 


10 ^ 0 . 30103 pt„ 5 ^r/Po 

0 


if Pt,s,r ^ Pb: 

Otherwise. 


(1.5) 


The exponent is chosen such that if pt,s,r/ Po = 1 the detection probability Pt^s,r = 0.5. This 
is important to be consistent with the definition of the range of the day. 

For diffuse sensor models we need a different approach to evaluate the total detection prob¬ 
ability. Thus, we define a more general version of Equation (1.3) to compute Pt regardless 
of the chosen sensor model assuming all probabilities are independent as 


^r=i- n (1.6) 

{s,r)eSxR 

In order to account for different weighted targets, we introduce Vt as the value of target 
t. Hence we denote the expected reward for target t as VfPt- If all targets have the same 
value, we simplify the expected reward for target t by using Pf. Based on this, we can 
define multiple objective functions depending on the chosen sensor model. While maxi¬ 
mizing the total expected reward is our only objective when using a cookie cutter model, 
we can also choose to maximize the minimum expected reward. The objective functions 
are summarized in Table 1.1. 


Objective 

Sensor Model 

Formula 

maximize total expected reward 

cookie cutter 

max Y, VtPt 
teT 

maximize average expected reward 

diffuse 

max Y VtPt/\T\ 

t^T 

maximize minimum expected reward 

diffuse 

max min V;/) 

teT 


Table 1.1: Objective functions - Three different objective functions are possible depending on 
the chosen sensor model, where V; denotes the value or weight of target t. 


1.2.2 Deployment and Usage 

There are different approaches to quantify the effectiveness of the deployment and usage 
of multistatic sensors. A strategic rather than a tactical approach is analyzed by Washburn 
(2010). He assumes that sources and receivers are deployed uniformly at random within 
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some region. The author develops methods for approximating the detection probabilities 
as a function of number of sensors. The most interesting aspect about this approach is the 
fact that it does not need to consider the geometric arrangement of the sensors. 

The DSTO analyzes multiple scenarios. Fewell and Ozols (2011) compare a field of mono¬ 
static sonar systems with that of a field of similar sonars operated multistatically, where 
sources and receivers are collocated. The direct comparison of the two modes of operation 
reveals that the correct choice of sensor models affects the outcome. Using a cookie cutter 
sensor model, the researchers find no advantage to the MSN. However, using the expo¬ 
nential model, the researchers achieved the same results in an MSN with about one quarter 
the number of sensors as the monostatic sensors. This finding is important since it reduces 
the number of pings in a given field. As soon as a ping is sent out, a hostile submarine 
knows the location of the source and will depart, which consequently makes it harder to 
detect. Another reason to reduce the number of pings is the artificial stress for sea dwellers 
produced by sonar systems. 

There are three main types of sensing coverage problems in the literature: area search, 
barrier search, and point coverage. Using the same sensor models as in (Fewell & Ozols, 
2011), Ozols and Fewell (2011) compare various layouts for sensor positions in a MSN, 
where a large area has to be covered. From the 27 layouts tested, four are recommended for 
use depending on the ratio of source and receiver cost. Additionally, in a classified report, 
Fewell, Ozols, and Rzetelski (2011) discuss the deployment of MSNs for a barrier search 
problem. Gong, Zhang, Cochran, and Xing (2013) address the barrier search problem by 
finding the optimal deployment of multistatic sensors to maximize the worst-case intrusion 
detectability. The authors show that a balanced structure yields the optimal solution. 

Craparo and Karata§ (2014) engage the point coverage problem for MSNs. The idea behind 
is to monitor a number of points of interest (Pols), e.g., oil platforms, aircraft carriers, 
ports, etc., in order to detect hostile underwater assets. The authors assume fixed Pols 
and receivers, and discuss various approaches to optimally place multiple sources. To our 
knowledge this is the only existing study considering the point coverage problem for MSNs. 
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1.3 Objectives 

From the three main types of sensing coverage problems (barrier search, area search and 
point coverage), this thesis focuses on the point coverage problem for MSNs. The overall 
question is: where do we place sources and receivers such that the Pols are optimally 
covered? To be in accordance with the terminology of other literature, ‘Pols’ are called 
‘targets’, and instead of ‘covering’ we say ‘detecting’. To engage this problem, we consider 
the following questions. 

• What are basic observations regarding point coverage sensing in MSNs? Particularly, 

- Is the convex hull spanned by the targets guaranteed to contain the optimal 
locations of sources and receivers? 

- Are there other ways to exclude suboptimal solutions? 

• How can we find an optimal placement in MSNs? Particularly, 

- Are there models or algorithms that represent point coverage sensing? 

- How do the developed means perform and how can they be improved? 


1.4 Scope, Limitations, and Assumptions 

For the developed models we assume that all targets and sensors are stationary and exist 
in a two-dimensional Euclidean plane with homogeneous environmental conditions. In the 
course of this study, all three detection probabilities shown in Figure 1.3 are considered. 
When we build our models, however, we start with the simple cookie cutter model and later 
advance to the diffuse sensor models. The effect of the blind zone can be greatly reduced 
by pulse compression as shown by Fewell and Ozols (2011, p. 11). For this reason, we do 
not consider blind zones in this thesis, i.e., we assume pb = 0. Furthermore, we neglect the 
target’s aspect dependence and assume all probabilities are independent. A set of notations 
used throughout the study is summarized in Table 1.2. 

1.5 Contributions And Outline 

In the first part of this thesis we explore some general observations about point coverage 
sensing with MSNs. We prove that a convex hull encasing the targets contains optimal 
sensor locations. Furthermore, we show the importance of various circles and the clusters 
of targets induced by their intersections for cookie cutter sensor models. We show that these 


7 



Element 

Set 

Deseription 

t 

T 

target 

s 

S 

souree 

r 

R 

reeeiver 

Vt 

M+ 

value of target t 

{xi,yi) 


eoordinates of objeet /, V/ G T U 5 U/? 


M+ 

Euelidean distanee between objeet i and j, \/i,j ^ TUSUR 

Po 

M+ 

range of the day 

Pt,s,r 

M+ 

equivalent monostatie range for target t, souree s, and reeeiver r 

Pt,s,r 

M+ 

probability to deteet target t with souree s and reeeiver r; Pt^s,r < 1 

P, 

M+ 

total deteetion probability for target t; P? < 1 


Table 1.2: Notation - This table summarizes the notation used throughout the study. 


clusters can be used to define bounds on the sensing problem. Additionally, we discuss how 
integer nonlinear programs (INLPs) can be used to find an optimal solution. 

The second part focuses on the development and enhancement of a new algorithm to find 
the optimal source position, assuming fixed targets and receivers in place. We investigate 
many details of the algorithm and assess methods and means to improve it. 
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CHAPTER 2: 

Observations On Point Coverage Sensing 


This chapter analyzes the point coverage scenario and discovers some basic observations. 
Those observations help to limit the solution space and exclude infeasible or non-optimal 
solutions. We also show how to apply the gained insights and whether the outcome justifies 
utilizing the developed method. 

2.1 The Convex Hull 

A simple approach to limit the set of potential sensor locations would be defining a rectan¬ 
gle that contains all Pols. Even though this can easily be implemented, we want to exclude 
as many locations as possible. A more sophisticated approach involves looking at the con¬ 
vex hull spanned by the targets. As vividly described in (De Berg, Van Kreveld, Overmars, 
& Schwarzkopf, 2000, p. 3), one can imagine the targets as nails sticking out of the plane. 
If we hold an elastic rubber band around the nails and let it go, the enclosed area will be 
the convex hull. 

2.1.1 Properties Of The Convex Hull 

The example in Figure 2.1 shows a convex hull for a particular set of targets. It also shows 
that all targets are either inside the hull or at one corner of the resulting polygon. The latter 
are called vertices. Hence we define the set C as 

C = {tGr|tisa vertex of the convex hull of T}. 

The convex hull is defined as the set Conv(r), that contains all points, that are spanned 
by the vertices in T. It is easy to see that Conv(r) = Conv(C). If all nails that are not 
vertices are removed in the previously mentioned example, the rubber band would still 
enclose the same area. Alongside this picturesque description of a convex hull, there also 
exists a mathematical definition. Each point p G Conv(C) with coordinates {xp,yp) can be 
written as a linear combination of the vertices in C. This is shown in Equations (2.1), where 
0 < Ar < 1 for all t G C. 
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Figure 2.1: Convex Hull Example - A Convex Hull spanned by a set of targets. Its vertices are 
Ti, T(i, Tg, Tio and Ti. 


X 


P ~ 


teC 


yp = 

leC 


1 = IA,. 

teC 


( 2 . 1 ) 


Theorem 2.1 deseribes the relationship between the eonvex hull and the deteetion proba¬ 
bility for sensor p, where Pt{p) denotes if p is a souree and Pt^s,p otherwise. 

Theorem 2.1. For every sensor position p ^ Conv(C), there exists a position p' G Conv(C), 
such that Pt{p) < Pt{p')yt G T. 


Proof. Assume sensor position p ^ Conv(C). Let p' G Conv(C) denote the position inside 
the eonvex hull with the shortest distanee to p. Then p' is either on an edge of the eonvex 
hull or eolloeated with one of its vertiees as shown in Figure 2.2. In both oases, the eonvex 
hull and with it all targets t G T are behind an imaginary line perpendioular to the line from 
p to p' orossing that line at p'. Otherwise there would be another p' that is eloser to p, 
whioh is a oontradiotion. 
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Figure 2.2: Theorem 2.1 Proof - There exist two cases for the closest position p' e Conv(C) to 
p. The left figure displays p' on an edge of the convex hull, while p' is collocated with a vertex 
on the right. In both cases all targets f e T are right of the dashed line, which is perpendicular 
to the line from p to p' and crosses p'. 


Let t be a target inside the convex hull. Without loss of generality, we assume that p and 
p'dSQ lying on a line parallel to the x-axis. Then the horizontal distance between p and t is 
\xp —Xt\ = \xp— Xp !I + \xpi — Xt I and we derive 

dip = {^p-^tf + {yp-ytf 

= i\xp-Xpi\ + \xpi-Xt\f+ {yp-ytf 
> [xp! -Xtf + {ypi -jtf 



Let pt{p) denote p?,p,r and pt^s,p if p is a source or receiver respectively. Thus the equivalent 
monostatic range for p' is 

Pt{p') = \ldt,p’dt^x < \/dt,pdt^x = Ptip), 

where x is the second type of sensor needed in an MSN. It follows that the detection 
probability for a target t using sensor p', Pt{p') < Pt{p) since all sensor model functions 
are monotonically nonincreasing with distance as is shown in Figure 1.3. □ 


We conclude that applying Theorem 2.1 we are able to place all sources and receivers 
inside the convex hull without sacrificing detection probabilities. Hence, we are now able 
to constrain our solution space, i.e., we can find optimal positions for each sensor inside 
the convex hull. We now describe and assess two ways to construct the convex hull. 
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2.1.2 A Linear Program To Find Vertices 


The first method uses a linear program (LP) to find the vertices of the convex hull. We use 
the fact that vertices themselves cannot be written as a linear combination of other points, 
whereas targets that are not vertices can. Based on this insight, we can define the following 
LP. 

Indices and Sets: 

t,t’&T targets. 


Data [units]: 

Xt X coordinate of target t [unitless], 

yt y coordinate of target t [unitless]. 

Decision Variables: 

A/ coefficient for target t's position subject to target t. 

Formulation: 


min 

A 


teT 


Subject to 



Xt' = ^ 
teT 

Vt' e r, 

(2.2a) 

yf = E ^t'yt 

teT 

'it' e T, 

(2.2b) 

1 = IA/' 

teT 

it' e r, 

(2.2c) 

0 < A/' 

itj' G T. 



The LP tries to minimize the number of target locations not expressed as a convex combi¬ 
nation of other target locations. For a non-vertex target t A/ can be zero, since t’s position 
can be represented by a linear combination of the vertices. The only possible value for 
vertices, however, is one. Thus, we can define the set of vertices C as 

C = {teT\Xl = \}. 
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Implementing and solving this LP in General Algebraic Modeling System (GAMS) runs 
very quickly for small numbers of targets. Figure 2.3, however, shows that more targets 
result in an disproportional increase in run times when using the simplex algorithm. In 
fact, Klee and Minty (1972) show that in a worst-case scenario, the simplex algorithm can 
reach exponential complexity. 



Figure 2.3: Convex Hull LP Run Time - For small numbers of targets the LP runs very quickly. 
Run Times, however, grow very quickly with the number of targets. 


2.1.3 The Graham Scan Algorithm 

Another algorithm to find the vertices of a convex hull with complexity 0(|r|log|r|) is 
the Graham scan described in (Graham, 1972). The algorithm starts with finding to G T, 
the target with the smallest value for yt. If there are multiple targets that share the smallest 
yt, we pick the one with the smallest Xt out of the candidates. 

Next we compute the angles Of each target t G T makes with to and the x-axis using the 
formula 


Ot = atan ^{yt -yt^.Xt-Xt^). (2.3) 

Sorting the targets by 9t starting with the smallest, the algorithm now considers each target 
as possible vertex. At each step it is determined whether the current target t and its two 
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predecessors t —I and t — 2 make a turn clockwise or counter clockwise by calculating 


C = {xt-i -Xt-2) X (yt-yt-i) - (yt-i -yt-2) x {xt-Xt-2)- (2.4) 

If the result is positive the points make a counter clockwise turn and the next target in the 
ordered list is considered. If c < 0 then the turn is counter clockwise and the middle target 
t — 1 is identified as non-vertex. At c = 0 all three targets are collinear, which also leads to 
discarding t — 1. Figure 2.4 displays the working of the Graham scan algorithm. 





Figure 2.4: Graham Scan - The working of the Graham scan is displayed. Targets are considered 
in order of the angle they make between the lowest target and the x-axis. In this example at the 
second step the algorithm detects a clockwise turn, which leads to discarding the middle target 
on step 3. 


The low complexity of the Graham scan is a result of the simple formula in Equation (2.4) 
to decide the complex issue about clockwise or counterclockwise turns. Outputting the 
vertices of the convex hull in a counter clockwise order is another advantage we later exploit 
in Section 3.4.3. 


2.2 Range Of The Day Circles 

This section analyzes the relationship between the range of the day, po and the cookie cutter 
sensor model. We already know from Equation (1.2) that a target is detected if and only 
if its equivalent monostatic range pt,s,r < Po- Besides that, we introduce range of the day 
circles (RDCs), i.e., circles around targets with radius po as shown in Eigure 2.5. In the 
following, we discuss methods that use RDCs that bound the problem. Even though we 
develop these methods for cookie cutter sensor models, they also apply to problems where 
we have to meet a particular detection probability, e.g., 80%. Here we simply change the 
radius to a value that corresponds with the demanded probaility. 


14 






Figure 2.5: Range Of The Day Circles Example - Assuming a cookie cutter sensor model, 
range of the day circles determine bounds for the problem. Each circle has a target for center 
and po for radius. 


2.2.1 Properties 

An important observation is the relationship between sensor positions relative to RDCs and 
a target’s detection probability. 

Theorem 2.2. A target t is detected only if at least one sensor is inside t ’s RDC. 


Proof. To arrive at a contradiction, we assume all sensors are located outside the RDC of 
target t and Pt^s,r = 1 for some 5 G 5 and r G R. Then dt^s > Po, V 5 G S and dt^r > Po, Vr G R. 
It follows that 


Pt,s,r — y/ dt^s X dt^r 

> a/Po X po = po- 

But then according to Equation (1.2) ^ ^ = 0, which is a contradiction. □ 

This, however, does not mean that we can find an optimal sensor placement by solely look- 
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Figure 2.6: Range Of The Day Circles Properties - The left plot illustrates that an optimal 
placement can be found by placing the source Si outside RDCs. Forcing the Si inside the circles 
results in losing coverage of targets. The right plot shows an example where an optimal receiver 
position is outside RDC intersections. Here, it is not possible to cover all targets by placing all 
receivers into RDC intersections. 


ing at positions inside RDCs. The left plot in Figure 2.6 demonstrates a eounterexample 
where the optimal position for the souree is loeated between both eliques of targets out¬ 
side any RDCs. Moving this souree inside a RDC of one of the eliques results in losing 
eoverage of targets in the other elique. 

A reasonable heuristie may involve ehoosing to plaee reeeivers into the eireles and let the 
sourees be the eonneeting elements of the MSN in order to keep the number of sourees 
as small as possible. However, the impression that optimal reeeiver positions ean only 
be found in RDC interseetions, unless a eirele does not have an interseetion with another 
eirele, is fallaeious as pietured by the right plot in Figure 2.6. Reeeiver Ri eannotbe moved 
into the interseetion of Ti’s and T 2 S RDCs without saerifioing eoverage of T\. Covering 
targets Tt, and inhibits the souree from moving eloser to T\. We take a eloser look at this 
example in Seetion 2.2.3. 


2.2.2 Clusters 

Nonetheless we ean use RDCs and their interseetions to form elusters. A eluster G C T 
is a maximal set of targets, where the RDCs of all targets t G G have at least one point in 
eommon. Sinee we require a eluster to be maximal, it is not possible to add another target 
to a eluster. In our example from Figure 2.5 we ean find the elusters 


G2 = {W2}, G3 = {T2,Te,Tg}, 

G5 = {T4,Ts}- 


Gi = {Ti}, 

G 4 = {T 4 ,Ts,Tiq}, 
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We make some observations about clusters. First of all, each target t G T occurs in at least 
one cluster. If there are no intersections with circles from other targets, a target forms its 
own cluster, e.g., Ti is the only element of G\. Additionally, it is possible that a target is a 
member of multiple clusters, e.g., T 4 G G 4 and Moreover, it is important to notice 

that even though at a first glance it looks like the clique of targets on top of Figure 2.5 (r 4 , 
T^, Tg and Tiq) all belong together, they actually form two separate clusters since Tg is only 
connected to r 4 . 

Figure 2.7 shows a special case: the left plot forms exactly one cluster, while the right plot 
forms three clusters {ri,r 2 }, {r 2 ,r 3 } and {ri,r 3 } since the three targets do not have a 
single point in common, even though the targets’ RDCs mutually intersect. 



Figure 2.7: Clusters Special Case - The targets in the left plot form exactly one cluster. The 
targets on the right, however, form the three clusters {ri,r 2 }, {T 2 ,T^} and since the 

three targets do not have a single point in common. 

Based on this, we define the minimal set of clusters G as the smallest set of clusters that 
contains all targets t G T. For the right plot in Figure 2.7 a minimal set of clusters is 
G = {{Ti, r 2 }, {r 2 , r 3 }}. Furthermore, we formulate the following theorem. 

Theorem 2.3. A lower bound on the number of sensors required to detect all targets is |G|. 

Proof. To arrive at a contradiction, we assume that all targets are covered with |G| — 1 
sensors. It follows from Theorem 2.2 that at least one sensor has to be inside each target’s 
RDC. Since the circles of each cluster have at least one point in common, putting a sensor 
on that point results in having a sensor in each RDC of this cluster. Moreover, since clusters 
are maximal, we cannot add another target. Hence, one cluster and therefore at least one 
target remains without a sensor in its respective RDC. □ 
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As a result from Theorem 2.3 we ean define |G| as a lower bound for the number of sensors 
required to deteet all targets. By the same token, we are able to find an upper bound for the 
number of targets we ean eover with a given number of sensors n by ehoosing n sets from 
G, sueh that the number of eontained targets is maximized. 

2.2.3 Other Circles 

Next to RDCs, there exist other eireles that provide insight into where to plaee either 
sourees or reeeivers. Craparo and Karata§ (2014) introduce the notion of the detection 
disk for each target t and receiver r that has t for center and its radius defined by 



In order to better distinguish detection disks from a second type of circles we introduce in 
this section, we refer to them as receiver circles (RCs). In order to detect a target t, we have 
to place a source inside fs RC. 

Figure 2.8 shows the example from Figure 2.6 augmented with its RCs. Here the red and 
blue circles are RCs for Ri and R 2 , respectively. Using Ri and one source we are able to 
detect Ti and T 2 simultaneously since their RCs overlap. If we want to detect Tt, or T 4 , 
however, we are not able to detect more than one target since their Ri RCs do not overlap 
with other Ri RCs. On the other hand, the R 2 RCs of and T 4 overlap and even have a 
common point with the RC of T 2 , such that placing a source here results in detecting T 2 , Tt, 
and T 4 with Ri. 

By the same token, we can define source circles (SCs) whose radius is determined by a 
source location as 


Figure 2.9 displays the SCs in our previous example. With the shown fixed source position, 
placing a receiver at the intersection of T^’s and r 4 ’s SCs ensures detecting both targets as 
well as T 2 . Placing a second receiver somewhere into Ti’s SC always detects T 2 also since 
T\s circle is completely overlapped by r 2 ’s SC. 
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Figure 2.8: Receiver Circles Example - Red and blue circles represent receivers i?i's and 7?2's 
receiver circles respectively. Placing a source at RC intersections ensures detection of respective 
targets. 



Figure 2.9: Source Circles Example - Placing a receiver at the intersection of targets T^'s and 
T 4 ’s source circles guarantees detecting both targets. A receiver somewhere in Ti's SC always 
detects T 2 as well because of the total overlapping. 


Assuming the positions of one type of sensors is already fixed, utilizing RCs and SCs is a 
way to decide positions for the remaining type. Nevertheless, a large number of targets and 
sensors makes their manual use hardly practical. 
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2.3 Utilizing Nonlinear Programs 

We now investigate the usage of INLPs to find optimal sensor placements in MSNs. The 
first model finds an optimal solution for a bistatic sensor network, i.e., exactly one source 
and one receiver, assuming a cookie cutter sensor model. Then we discuss how INLPs 
perform on extended networks. 


2.3.1 Bistatic Cookie Cutter Sensor Networks 

For this model we start from the premise that we already found the set of vertices C, using 
the Graham scan, for example. Then we formulate the following INLP to find optimal 
positions for a bistatic sensor network assuming a cookie cutter sensor model. 

Indices and Sets: 

t G r targets, 

p & P = {s,r} type of sensor, 5 for source and r for receiver. 

Data [units]: 

Xt X coordinate of target t [unitless], 

yt y coordinate of target t [unitless], 

Vt weight of target t [unitless], 
po range of the day [miles], 

C vertex set, C = {t G r|t is a vertex} [unitless], 

Mt penalty for not covering target t [unitless]. 

Decision Variables: 

Xj’ coefficient for sensor p’s position subject to vertex t, 

Xp X coordinate of sensor p, 
yp y coordinate of sensor p, 
ht binary, 1 if target t is not detected, 0 otherwise. 
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Formulation: 


min Y^Vfht 
Subject to 

leC 

jp = E yt 

teC 

1 = E V 

teC 

po "b Mihi ^ dt^s X dpf 

0<xf 

hte {0,1} 


Vjo G P, 

'ip G P, 

ip G P, 

Vt G T, 
it G G P, 
Vt G P. 


(2.5a) 

(2.5b) 

(2.5c) 

(2.5d) 


The model forces the positions of the source and receiver to be inside the convex hull 
by defining them as linear combinations of the vertices in Equations (2.5a) to (2.5c). We 
introduce the binary decision variable ht for each target t G P specifying in Equation (2.5d) 
whether a target t is detected, ht = 0, or hidden, ht = 1. If t is not covered, then a large 
number Mt is added to the range of the day, po, making sure the constraint is still feasible. 

Camm, Raturi, and Tsubakitani (1990) show that tight penalty parameters like Mt reduce 
the time a solver needs to find an optimal solution. In order to make sure the INEP is always 
feasible, the left-hand side (EHS) of this equation has to be equal to the largest product of 
distances possible for each target. Since the largest distance from a target position to a 
point in a convex hull is to one of its vertices, we can set Mt to its square. We can also 
substract po since it is already added on the EHS. Thus we have 


Mt = rmxd}fi 

t’ec ’ 


-Po- 


( 2 . 6 ) 
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It is important to notice that there is generally more than one optimal solution with this 
model. First of all, interehanging the souree and reeeiver in a bistatie environment does 
not ehange the outeome. The Cassini oval drawn by both sensors is still the same. But 
more importantly, if a sensor is moved slightly in one direetion, it may still deteet the same 
targets. And last but not least, there might be a solution deteeting different targets, but the 
same number of targets. Henee, in most seenarios there exist infinitely many solutions. 

2.3.2 Other Networks 

The model deseribed in Seetion 2.3.1 only supports bistatie eookie eutter sensor networks. 
Consequently, there are two ways to extend this model: finding positions for multiple 
sourees and reeeivers, and expanding to diffuse sensor models. Though models ean be 
formulated, up to this point we are not able to define a model that runs in reasonable time. 
In a future study, advaneed versions of this model ean be further investigated. 
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CHAPTER 3: 

Divide Best Sector Algorithm 


In this chapter we slightly change the scenario. In an operational setting, it is not unlikely 
that receivers are already deployed throughout a field by plane or ship. As soon as a hostile 
underwater asset is suspected in the area of responsibility (AoR), e.g., by passive detection, 
sources are used to send out pings in order to find the exact position of the intruder (Navy 
personnel, personal communication, April 2014). As discussed, it is crucial to send out as 
few sources and subsequent pings as possible. 

3.1 Problem Statement 

Consider the following scenario: there is a set T of targets with fixed positions that have 
to be monitored or detected. Furthermore, there is a set R of receivers whose positions we 
already fixed. The objective is to find optimal positions for a limited number of sources, 
defined in set S. We do not assume a particular sensor model; rather, we discuss a method 
that works with all sensor models shown in Figure 1.3. 

Figure 3.1 shows a generic scenario with fixed targets and receivers. The z-axis represents 
the average detection probability assuming a single source is deployed at the respective 
position. Multiple local maxima in this plot show that the objective function is nonconvex. 
Nonetheless we now develop an algorithm guaranteed to solve this problem to near global 
optimality when |5| = 1. 

3.2 Algorithm Development 

In this section we develop the Divide Best Sector (DiBS) algorithm. The algorithm par¬ 
titions the area of possible solutions into sectors. At each iteration the sector with the 
highest upper bound for the objective function is further divided into smaller sectors. In 
Section 2.1 we showed that optimal sensor positions can always be found inside the convex 
hull spanned by the targets. We slightly relax this condition to the smallest rectangle with 
edges parallel to the .r- and y-axes that contains all targets. This way we can easily divide 
the area of possible source locations into smaller sectors by slicing horizontally as well as 
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Figure 3.1: Problem Statement - The optimal source position with fixed target and receiver 
locations needs to be found. The z-axis represents the average detection probability assuming a 
source is deployed at the respective position. Multiple local maxima indicate nonconvexity. 

vertically. 

Each sector 7 of the resulting set of sectors E now is evaluated with respect to an up¬ 
per bound, uB( 7 ), for the objective function. This is the heart of DiBS and illustrated in 
Figure 3.2. Since receivers and targets occupy fixed positions in this scenario, the target- 
receiver distance dt^r is constant. Thus, the equivalent range pt,s,r and subsequently the 
detection probability Pt only depend on the target-source distance dt^s- Hence, for every 
target we determine the hypothetical source position that is closest to the target but still 
inside 7 . These hypothetical source positions are either on the edge of the sector or the 
target position itself, as seen in Figure 3.2. 

With the shortest distance dt^j from target t to its respective hypothetical source position 
inside sector 7 we can compute the highest possible detection probability Pt, assuming the 
source location is inside 7 . Having done this with all targets, we compute the upper bound, 
uB ( 7 ), according to Table 1.1. The following pseudocode finds the upper bound for a given 
sector 7 , where x™" (y™'") is the lowest and the highest x (y) coordinate for 7 . 
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Figure 3.2: Determine Upper Bound - For each target t there exists a hypothetical source 
position inside sector / that has the shortest distance dtj to the target. This distance determines 
the highest possible detection probability with a source position inside y. Merging the probabilities 
for all targets results in an upper bound for y. 


1: procedure Upper.Bound( 7 ) 

2 : for allt eT do 

3: ^ min (max (x™", Xf), 

4: ^ min(max(3;™'',yf),y™“) 

5: dt^Y ^ \/{^t-XsY + {yt-ysY 

6 : compute Pt 

7: end for 

8: compute objective value Z 

9: return Z 

10 : end procedure 


> use chosen sensor model 
> use chosen objective 


In a next step we pick the sector 7 G F in which uB( 7 ) is maximized and divide it into 
smaller sectors with respective upper bounds. Consistently repeating this method results in 
smaller sectors and subsequently tighter upper bounds until a termination condition is met. 
Various termination conditions are conceivable, such as a maximum sector size, optimality 
range, etc. The following pseudocode shows the workings of DiBS. 
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1; procedure DiBS 
2 : create initial F 

3: optimal ^ FALSE 

4: while -^optimal do 

5; select 7 G r such that Upper.Bound(7) is maximal 

6 : if termination condition is met then > use chosen termination condition 

7; optimal ^ TRUE 

8 : else 

9: r' ^ split 7 into smaller sectors 

10: r^rur\{7} 

11: end if 

12: end while 

13: return 7 

14: end procedure 

Eigure 3.3 shows the example from Eigure 3.1 solved with DiBS. It vividly illustrates how 
sectors with a low upper bound are ignored while promising areas are further explored. Eor 
the rest of the chapter we look at various details of DiBS and give recommendations on 
how to use it most efficiently. 


3.3 Termination Conditions 

Eirst, we want to make sure, is that the algorithm eventually finishes when it meets a specific 
condition. This section discusses various termination conditions based on sector size and 
optimality gap. 

3.3.1 Termination By Sector Size 

There are two ways to define the sector size: by its area and by its edge lengths. A small 
area, however, can be the result of a very narrow but long rectangle. This metric is not very 
useful, because detection probabilities can change greatly along a line, as well as in a long 
rectangle. Thus, even if a sector with small area has a high upper bound, a sensor placed in 
this sector may perform badly. Additionally, telling a decision maker to place the source in 
a 0.01 X 10,000 foot rectangle might not be useful in practice. 
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Figure 3.3: Algorithm Example - A scenario with ten targets and three receivers is given (see 
Figure 3.1). DiBS finds the source position with the highest average detection probability by 
dividing only promising sectors, i.e., sectors with a high upper bound, into smaller ones. 


On the other hand, using a maximum edge length strongly eonstrains possible souree loea- 
tions. Choosing a maximum edge length of ten feet, for example, results in a final seetor 
that has size of at most 10x10 feet, a suffieient preeision for helieopter dips. Moreover the 
ehange in deteetion probabilities in sueh a region is very limited. 

Therefore, we disregard using area as the termination eondition. For the rest of this ehapter, 
termination by seetor size denotes termination when the length of the longest edge of the 
seetor with the highest upper bound is suffieiently small. 

3.3.2 Termination By Optimality Gap 

Though simple and useful, it turns out that termination by seetor size has a drawbaek. If 
we request a seetor with very short edges, the differenee for the upper bounds between 
adjaeent, same-sized seetors beeomes negligible. This leads to dividing all seetors with the 
same size in an area before eontinuing with smaller ones. This, however, takes the majority 
of the time and, on top of this, the gain in objeetive value is almost always negligible. 

Henee, we expand the termination eondition with a teehnique used by solvers: we introduee 
a tolerated optimality gap, £. By allowing DiBS to stop as soon as it finds a speeifie sensor 
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position s' that has an objective value within e of the upper bound, we are able to stem 
the described effect. The effectiveness of the optimality gap strongly depends on which 
position we choose for s'. 

We take a closer look at the methods summarized by Table 3.1. A method with fewer 
potential source positions, like Center, does not require high computational power but is 
subject to missing the optimality gap more often. On the other hand, checking five positions 
with Corners+Center is more precise but increases run time. The last method. Targets, is 
an add-on for the previous ones. The idea is that the detection probability for a single target 
is replacedmaximizedmaximum only if at least one sensor is at the target’s position. Since 
there are sectors that do not contain any targets, it can only be used in conjunction with 
another method. 


Method 

Amount 

Remarks 

Center 

1 

uses seetor’s eenter point 

Corners 

4 

uses seetor’s vertiees 

Corners-i-Center 

5 

eombines Corners and Center 

Targets 

0,1,... 

uses target positions inside seetor 


Table 3.1: Optimality Gap Methods - These Methods look for a specific source position inside 
a sector in order to determine the optimality gap. The last method, Targets, has to be used in 
conjunction with one of the other methods since sectors with no targets inside exists. 


The example in Figure 3.4 shows how an optimality gap limits the number of created and 
evaluated sectors and thus the run time. In this particular problem instance the baseline 
case created more than 16,000 sectors for a demanded precision of 0.001. Utilizing an 
optimality gap of £ = 0.05, however, already terminates DiBS after about 200 created 
sectors with only a minor difference between Center and Comers. 

Running more examples, we observe that the gap between Center and Comers always is 
small. Most of the time Comers terminates with fewer sectors than Center, because it can 
choose from four potential source positions instead of just one. Center, by contrast, has 
the better run time. Since Corners-i-Center combines both methods it always uses the same 
number of sectors as the method with fewer sectors, but it takes the longest run time. In 
all tests we ran, we never observed that adding the method Targets had any effect on the 
outcome except for increasing run time. 
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Figure 3.4: Optimality Gap Example - While the number of evaluated sectors for the baseline 
case without optimality gap feature increases with higher demanded precision, using an optimality 
gap terminates the algorithm earlier. 

Including the termination eondition with an optimality gap is strongly reeommended sinee 
it signifieantly reduees the number of created and evaluated sectors, especially for high 
precision calculations. Even though the run time per iteration slightly inereases, the total 
run time ean reduee to a fraetion of the original run time, whereas the ehoiee of method for 
finding the definite souree position s' hardly ehanges the outeome. 

3.3.3 Utilizing Lower Bounds 

On eloser examination, we realize that the methods described in the previous section pro¬ 
vide a lower bound, IB( 7 ). This section discusses whether we ean use lower bounds to 
disregard seetors. Additionally, we show the performanee of DiBS using lower bounds. 

Consider disearding a hypothetieal seetor / whose upper bound, uB(/), is less or equal 
than the lower bound, IB( 7 ), of another seetor 7 . Surprisingly, this is not neeessary. Sinee 
IB ( 7 ) is the objective value for a definite souree position in 7 , it is always included in one 
seetor of 7 ’s split. Henee there is always a seetor out of 7 ’s split that still has a higher upper 
and lower bound than /. Thus, / will never be chosen as the best seetor, and therefore, 
disearding does not reduee the total number of sectors created. 

It does, however, reduce the number of sectors in a list that has to be ordered by upper 
bounds. Most programming languages use fast sorting algorithms like Quicksort. So or- 
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dering a large list eompared to other operations in DiBS is negligible. 

Nevertheless, we ean use the lower bound to indireetly measure the goodness of the upper 
bound. A tight upper bound is important for faster identifieation of a solution. Figure 3.5 
illustrates the change of gap between upper and lower bound throughout the execution of 
DiBS. Here, the right plot shows the underlying trend of a decreasing gap the longer the 
algorithm runs. Jumps to lower values for the lower bound occur when the next sector 
with highest upper bound is much larger than the previous sector and therefore can provide 
a bigger gap between lower and upper bound. The left plot verifies our suspicion that 
smaller sectors yield smaller gaps and therefore tighter upper bounds. 
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Figure 3.5: DiBS Goodness - With increasing number of iterations the gap between upper and 
lower bound decreases (left plot). Steep cuts on lower bound values occur on iterations where a 
small sector is sufficiently split and DiBS jumps back to a larger sector and thus larger gap. The 
right plot mirrors the relationship between longest edge size as a multiple of longest initial edge 
and gap size. Smaller sectors yield tighter upper bounds. 


Summarizing, we conclude that lower bounds do not accelerate DiBS. They can, however, 
be used to show that upper bounds are getting tighter along the algorithm. Hence, the 
algorithm quickly achieves a small optimality gap, which can be used to terminate DiBS 
early as shown in Section 3.3.2. 


3.4 Sector Splitting 

Since the algorithm begins with creating the initial set of sectors F, this procedure and 
subsequently the dividing into smaller sectors is the next feature of DiBS into which we 
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take a closer look. The question is how many vertical and horizontal slices are efficient 
with respect to complexity and run time? First we investigate a fixed number of slices and 
expand later on to dynamic splitting methods. 

3.4.1 Fixed Splitting 

More slices reduce sector sizes faster and therefore meet the termination condition in fewer 
iterations of DiBS. Each iteration, however, requires more computations and thus more 
time. The plots in Figure 3.6 show the algorithm’s performance for a particular problem 
instance using various settings for initial division and further splitting. The derived source 
position is always the same for all settings. As assumed, the number of iterations decreases 
quickly with more initial sectors and even more with a higher number of new sectors per 
iteration. The positive effect, however, stagnates very quickly. 

number of initial sectors:-4; - 36; 100; 225; - 400 




new sectors per iteration new sectors per iteration 

Figure 3.6: Sector Splitting - Algorithm’s performance with different splitting settings. Splitting 
into many sectors per iteration increases the total number of sectors and subsequently the run 
time of the algorithm. The size of the initial set of sectors also affects run time especially for 
high numbers of new sectors per iteration. The reduction in number of total iterations is not as 
decisive. 

On the other hand, the right plot in Figure 3.6 reflects the increasing number of created 
sectors. While the change in the number of iterations is between 8 and 62, the values for 
the number of created sectors range from 270 to over 10,000. Running DiBS on a test 
machine shows that the run time is directly related to the number of created sectors, e.g., 
0.9 seconds for an initial and further split of 4 sectors versus 28.7 seconds for an initial and 
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further split of 400 sectors. The number of iterations does not seem to have any effect on 
the run time. 

The results are decisive for newly created sectors: dividing the best sector into fewer new 
sectors greatly decreases run time and improves algorithm efficiency. A reason for this 
observation lies in dividing small sectors. We know from Section 3.3.3 that small sectors 
have tighter upper bounds. Hence, the variance in upper bounds of their split sectors is 
small. So after splitting a small sector into a set of many small sectors F', it is very likely 
that the algorithm has to evaluate and, even worse, split most of the sectors in F' before it 
is able to continue with a smaller sized sector. This behavior creates many sectors that are 
unnecessarily evaluated and therefore waste resources. A different approach is discussed 
in Section 3.4.2. 

The creation of the initial set of sectors also has an effect on the total number of created 
sectors. Especially at high numbers of new sectors per iteration it becomes a critical factor, 
as seen in Figure 3.6. It appears a smaller initial set performs more efficiently. 

In most cases, starting with a small set of sectors and slicing the best sector into a small 
number of new sectors is most efficient and performs best. We recommend to always start 
with this setting and expand as circumstances require. 

3.4.2 Dynamic Splitting 

We discovered in Section 3.4.1 that splitting into many sectors per iteration produces a lot 
of overhead, especially at small sector sizes. The question arises whether one can reduce 
the number of created sectors via dynamic splitting. In dynamic splitting, large sectors or 
sectors with a high optimality gap are split into more sectors than small sectors or sectors 
with tighter upper bounds. We leave that question open for future work. 

3.4.3 Plane Rotation 

The DiBS algorithm starts with finding the smallest rectangle with edges parallel to the 
X- and y-axes that contains all targets. We know from Section 2.1 that the optimal sensor 
position lies inside the convex hull spanned by the targets. Hence, it is possible that the 
resulting rectangle contains a large overhang which is subject to being evaluated unneces¬ 
sarily. This section discusses the use of the minimum-area rectangle without constrained 


32 



edges and how to rotate the plane, such that DiBS still can be applied. 

Freeman and Shapira (1975) prove that the minumum-area rectangle encasing a convex hull 
has a side collinear with one of the edges of the convex hull, like the examples in Figure 3.7. 
Hence we choose the rectangle with smallest area out of the |C|, i.e., the number of vertices 
of the convex hull, rectangles we can form. Using the Graham scan to find C has the 
advantage that the vertices are already ordered counter clockwise. Therefore, we can use 
this output to easily construct the edges of the convex hull. 



Figure 3.7: Minimum-Area Rectangle Example - The plot shows the original rectangle as 
well as two rectangles that have one side collinear with one side of the convex hull. The Tg —Tg 
rectangle is the minimum-area rectangle that encases all targets. 

The example in Figure 3.7 shows the original rectangle produced by DiBS as well as the 
rectangles with the smallest and the largest area that have one side collinear with one of the 
edges of the convex hull. The respective areas in square units are 


original: 636, T(, - T^: 505.375, Tio - Ti: 619.420. 


This example shows that using the minimum-area rectangle can greatly reduce the over¬ 
hang. Having determined the minimum-area rectangle, we need to rotate the plane such 
that the edges of the rectangle are parallel to x- and y-axis. An efficient way to do this is 
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applying a linear transformation to the sets of target and reeeiver locations (Leon, 2010, 
p. 166ff.). Let to — h be the edge that is collinear with the minimum-area rectangle. Then 
the linear transformation is represented by the matrix 


A = 


cosO 
— sin0 


sin0 

COS0 


where 9 = tan ^ 


yti 

Xt, -XtJ 


The matrix multiplication AX provides the new coordinates for targets and receivers, where 
the first and second row of X holds the x and y coordinates respectively. 


Since plane rotation can have a huge impact on the size of the overhang, it should always 
be considered before running DiBS. There is also the fact that Graham scan and linear 
transformation are very efficient and therefore executed quickly. 


3.4.4 Handling Uneven Rectangles 

The next improvement on sector splitting aims for rectangles with uneven edge lengths. 
Obviously the short edge reaches its termination condition faster than the long edge if DiBS 
performs the same number of slices in both directions at each iteration. Since termination 
by sector size considers only the longest edge, the algorithm continues dividing sectors in 
both directions. Hence it creates more sectors than needed. 

There are two ways to address this problem: stopping splits in directions that meet the ter¬ 
mination condition and finding an initial split that has only square sectors. The first method 
is easy to implement by setting the number of slices for edges that meet the termination 
condition to zero. There is, however, already a problem with uneven rectangles before the 
short edge meets the termination condition. As opposed to squares with the same area, 
especially long rectangles are more likely to have a big variance in detection probabilities 
(see Section 3.3.1). This subsequently leads to looser upper bounds. 

The second approach bypasses this issue by starting with squares at the first iteration. We 
already discovered in Section 3.4.1, that a small set of initial sectors performs more effi¬ 
ciently. So we are looking for a method that creates a low number of squares that covers an 
uneven rectangle. 

Figure 3.8 shows two methods to create initial squares. For the rest of this section, we 
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Method 1 


% 


/-^ 





- J 


> 


Method 2 


a 


_/V_ 








b 


Figure 3.8: Creating initial squares - The first method uses half of the long edge b for square 
edge length, while method 2 uses the short edge. The red shaded area reflects the created 
overhang when using only squares for the initial split. 


assume that a and b are the short and the long edge respeetively. The first method outs 
the long side b in half and takes the result as the edge length for the squares. The seoond 
method, on the other hand, uses the short edge a as the edge length for the squares. Both 
methods, however, oreate an overhang, shaded red in the figure, that hold non-optimal 
solutions. The overhangs are oaloulated as follows: 

For method 1: — aZ?. (3.1) 

fO if a I h, 

For method 2: 0.2= \ (3.2) 

[[a — {b mod a)]a otherwise. 

In order to keep the unneeded area minimal, we always ohoose the method with the smallest 
overhang. Table 3.2 lists two examples with different ehoiees of the method used. 


a 

b 

r = a/b 

Method 1 Overhang 

Method 2 Overhang 

Select 

4 

5 

0.8 

52-4x5 = 5 

[4 - (5 mod 4)]4 = 12 

Method 1 

2 

5 

0.4 

52-2x5 = 15 

[2 - (5 mod 2)]2 = 2 

Method 2 


Table 3.2: Examples for Overhang - For the first example, method 1 has the smaller overhang, 
hence we select it. For the second example, however, we choose method 2 because of its smaller 
overhang. 


Going a step farther, we want to see whether we are able to decide which method to use by 
solely looking at the ratio between long and short edge. 
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Theorem 3.1. Method 1 has less overhang than method 2 if and only if ^<r <\, where 
r = alb. 


Proof Without loss of generality, we assume b =1. Then 0 < r = a < 1. Furthermore the 
overhangs for both methods compute as follows: 

Cli = b^ — ab = I — r, 

fO if a\b, 

^ 2 ={ 

[[a — {b mod a)]a = [r — (1 mod r)]r otherwise. 

Method 1 has less overhang if and only if < 0, which is visualized in Figure 3.9. 

With a ratio of r = 1, i.e., both sides have the same length, both methods do not have an 
overhang since a and b already form a square. Discontinuities occur on positions where a 
divides b, such as Yi, '/s, !4 and so on. Here, method 2 produces no overhang. Moreover, by 
inspection, we see — ^^2 > 0 for 0 < r < 0.5. 



Figure 3.9: Theorem 3.1 Proof - The figure shows the difference between method I's overhang, 
ill, and method 2's overhang, 0.2, with respect to the ratio r = alb. A positive value indicates 
a smaller overhang for method 2, while method 1 is preferred at negative values. 
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Hence, we are looking for a root in 0.5 < r < 1, where a\b and 1 mod r = 1 — r. Thus 


— ^^2 = 1 — r — [r — (1 mod r)]r 
= l-r-[r-(l-r)]r 
= l-2r2, 

of which the only nonnegative root is 77^- □ 


When DiBS splits a sector, it does not recognize whether or not a newly created sector is 
in the overhang. An additional check removes those sectors and prevents selecting them as 
best sector for further splitting. 


3.5 Finding Multiple Sensor Locations 

A way to increase the number of possible applications of DiBS is increasing the number 
of sensor positions optimized. This section investigates an iterative method as well as a 
simultaneous approach. 


3.5.1 Placing Sources Iteratively 

We consider a set S of sensors, that has to be optimally deployed in a scenario with fixed 
targets and receivers. Furthermore, we assume a cookie cutter sensor model. The position 
for the first source is determined by run DiBS without modifications. In a next step we 
remove all targets that are detected by the algorithm’s solution and running DiBS again for 
the second sensor position. We repeat those step until either all targets are detected or all 
sources are set. The following pseudocode illustrates the iterative method. 
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1; procedure Iterative.Method 
2: while |r| > 0 A |5| > 0 do 

3: select s e S 

4: run DiBS 

5; T ■^T\ {detected targets} 

6 : 5 ^ 5\{4 

7: end while 

8 ; end procedure 

Since we place sources in a greedy approach, the resulting solution generally is not optimal. 
Notwithstanding, this algorithm reveals some information about the optimal solution. First 
of all, in case all targets are detected by this method, the solution indeed is optimal. 

Otherwise, we can use the outcome to define a lower and upper bound for this problem 
by applying an observation made by Craparo, How, and Modiano (2011). Leveraging sub¬ 
modularity in a similar coverage setting, the authors determine that an iterative algorithm 
has an optimality gap of at most 1 — e^^ ~ 0.632. Hence, we define the lower bound, IB, 
as the outcome of the iterative approach and the upper bound 


^ [1.582 xlBj. 

We conclude that even though the iterative method does not guarantee an optimal solution, 
it is a quick way to get an estimate about the quality of other solutions. Combined with our 
observations from Section 2.2, we are able to narrow down the optimal solution even more. 

3.5.2 Placing Sources Simultaneously 

In contrast to the previous method, a different approach is placing multiple sources simul¬ 
taneously. In doing so, we place each source in an arbitrary sector and evaluate upper 
bounds of occupied sectors. Merging those upper bounds yields the value for the objective 
function. In order to exhaust all possible solutions, we evaluate all combinations of sec¬ 
tors. We have to take into account that multiple sources might occupy the same sector and 
that the order of sources is indifferent. Hence, by definition the number of combinations is 
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calculated by combinations of multisets (Brualdi, 2010, p. 52ff.) as 



(bi + in-i)! 

V |5| J 

|X|!(|r|-l)! 


where |5| is the number of sensors and |r| the number of sectors in a particular iteration. 
The base ease with one souree represents the original DiBS where the number of eombina- 
tions equals the number of seetors, i.e., plaeing the souree onee in eaeh seetor. With two 
sources and four sectors we already have the following ten eombinations. 

l:{7i,7i} 2:{7i,r2} 3:{7i,73} 4:{7i,74} 5:{72,72} 

6:{72,73} 7:{72,74} 8:{73,73} 9:{73,74} 10:{74,74} 

The number of combinations grows quickly with the number of sectors and sources whieh 
is visualized in Figure 3.10. While plaeing just a few sourees simultaneously might still 
be manageable, using this approach with many sensors quiekly becomes an issue. This is 
aggravated by the fact that selecting and dividing multiple seetors per iteration inereases 
the number of seetors and subsequently the number of eombinations even more quiekly. 
A future study might find ways to reduce that number or find another method to deal with 
multiple sourees. 


number of sources:-1; - 2; - 3; - 5; - 10 



Figure 3.10: Sector Combinations - This plot displays the number of combinations that have 
to be evaluated with respect to the number of sectors and sources. The base case with one 
source represents the original DiBS. 
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CHAPTER 4: 
Conclusion 


We started by developing some insights that help bound the problem. We showed that plac¬ 
ing sensors inside the convex hull that encases all targets results in detection probabilities 
as least as good as outside the hull. This important observation limits the area where we 
have to search for optimal sensor positions regardless of the sensor model we use. 

Assuming a cookie cutter sensor model, we can utilize clusters of RDCs. The number 
of elements in the minimal set of clusters, G, represents a lower bound on the number of 
sensors required to cover all targets. In addition, with a fixed number of sensors, n, we 
are able to determine an upper bound on the number of targets by selecting n clusters from 
G, such that the number of contained targets is maximized. Other circles, such as RCs and 
SCs, are means to narrow down possible sensor positions assuming the location of one kind 
of sensor is fixed. 

In a subsequent step, we formulated an INLP to find optimal positions for a bistatic cookie 
cutter sensor network. Up to this point, expanding this model to multiple sensors or other 
sensor models does not provide a solution in reasonable time. Hence, we developed the 
DiBS algorithm. 

DiBS assumes that target and receiver positions are fixed and determines the optimal posi¬ 
tion for a single source regardless of the applied sensor model. By dividing sectors with the 
highest upper bound for the objective function, it eventually reaches a termination condition 
like maximum size of a sector’s longest edge or a maximum optimality gap. 

Further investigations of the DiBS algorithm’s details resulted in the following observations 
and recommendations. The number of new sectors created at each iteration as well as the 
number of initial sectors should be kept small in order to reduce the total number of created 
and evaluated sectors and subsequently the run time of DiBS. 

Rotating the plane such that the edges of the minimum-area rectangle encasing all targets 
are parallel to the x- and y-axis results in a smaller needless overhang. Because of the high 
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efficiency to find the minimum-area rectangle and to rotate the plane, this improvement 
should always be considered. 

Furthermore, we discussed two ways to handle uneven rectangles that can lead to creating 
and evaluating many unnecessary sectors. Both methods approach this issue differently. 
Their efficiency depends on the given scenario. We recommend, however, to apply one of 
the methods to reduce the algorithm’s run time. 

Finally, we illustrate two methods to apply DiBS to problems where we have to find mul¬ 
tiple source locations. The faster iterative method assuming a cookie cutter sensor model 
provides lower and upper bounds on the number of detected targets. The simultaneous ap¬ 
proach, however, quickly creates a large number of combinations that have to evaluated. A 
future study can develop ways to reduce the number of combinations or find other means 
to find multiple source locations using DiBS. 
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APPENDIX A: 

Excel Analysis And Decision Tool 


For a better understanding and visualization of the observations from Chapter 2 , we devel¬ 
oped an analysis and decision tool in Microsoft Excel. It provides a number of features for 
a bistatic cookie cutter sensor model, e.g., showing the Cassini oval. 


The Data worksheet shown in Figure A. 1 contains the scenario data such as target positions 
and range of the day. Furthermore, it shows the positions of source and receiver and deter¬ 
mines the detection probability for each target based on target-source and target-receiver 
distance. 



A 

B 

C 

D 

E F G H 1 J K L 

1 






2 


Ranee of the Day 



Otoject 

X 

V 

Src Dist 

Rev Dist 

detected 


3 


Po 

3.5 



T1 

-15 

-4.75 

18.6269294 

15.3121684 

FALSE 


4 






T2 

6.25 

-9.5 

13.9950884 

19.054199 

FALSE 


5 


Source 

X 

V 


T3 

-9.25 

-11 

18.2705364 

17.9669836 

FALSE 


6 


Sourcel 

1.6 

3.7 


T4 

0.75 

8.75 

5.12103505 

5.70306935 

FALSE 


7 






T5 

-4 

5.5 

5.88217647 

0.89442719 

TRUE 


8 


Receiver 

X 

V 


T6 

9 

-11.25 

16.681202 

22.0808175 

FALSE 


9 


Receiver 1 

-4.4 

6.3 


T7 

-9.5 

-15 

21.746264 

21.9020547 

FALSE 


10 






T8 

2.75 

11.5 

7.88431988 

8.84095583 

FALSE 


11 


UpdateTarget Information 

1 

T9 

3.5 

-12.25 

16.0627675 

20.1621551 

FALSE 


12 


J 

TIO 

-4 

7.5 

6.76756973 

1.26491106 

TRUE 


13 








14 


LoadTargetstrom .CSV 

1 








15 
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Figure A.l: MSN Tool Data Sheet - The Data worksheet contains information about targets, 
sensors and range of the day. It provides methods to change the scenario. 


There are two ways to change the target data: by changing the values in the table or by 
loading a comma-separated values (CSV) file. Using the first method, we have to press 
the Update Target Information button after all changes were made in order to apply the 
changes. The second method initiated by clicking Load Targets from .CSV does not require 
this step. The following example shows the required structure for a CSV file. 


11,-15,-4.75 

12,6.25,-9.5 
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Figure A.2 shows the Chart worksheet that visualizes the seenario. Here, we are able to 
show or hide various features. The serollbars on the top and on the right move either source 
or receiver depending on what is chosen on the left. Detected targets are represented by 
empty circles; otherwise, targets are represented by filled circles. The Data box contains 
some general information while the Selection box provides a short report about each ele¬ 
ment in the chart based on the mouse position. 


■ Options- 

[7 Show Cassini Oval 
W Show Convex Hull 
W Show Range of the Day Circles 
r* Show Soufce Circles 
r” Show Receiver Circles 
Use Scrollbars to 
O move Source ♦ 

(S move Receiver A 


• Data- 

Range of the Day: 3.5 
Source position: (1-6; 3.7) 
Receiver position: (-4.4: 6.3) 
No. of covered Targets: 2 


- Selection- 

— Range of the Day (T10) 

This circle intersects with the 
Range of the Day circle of: 

T4.T5 

In order to detect T10, at least one 
sensor has to be inside this circle. 


Find Optimum 



Figure A.2: MSN Tool Chart - The scrollbars move either source or receiver based on the 
selection on the right. The Selection box provides additional information about shown objects in 
the chart when hovering with the mouse. In this example target Tio’s RDC is selected. 

The Find Optimum button is without functionality in this version of the tool. Here, the 
INLP from Section 2.3 can be integrated in future work. 
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APPENDIX B: 

R Implementation Of DIBS 


DiBS <— function (targets , receivers , rho_0 , longest_edge , verbose = TRUE, obj = ’ avg ’ , 
sensor = ' fermi ’ , b = 0.25, gap_method = ’none’, gap = 0.05) { 

# Main function to run DiBS 

# 

i krgs : 

M targets: Data frame with object, x, y and value. 

^ receivers: Data frame with object, x and y. 

M rho_0: Range of the day. 
i longest_ edge: Termination condition. 

M verbose: If TRUE, prints progress. Default is TRUE. 

# obj: Objective. Possible values: *avg’, ^min\ Default is ^avg\ 

i sensor: Sensor model. Possible values: ^cookie*, *fermi*, ^exp ^. Default is ^fermi^. 
M b: Diffusivity parameter for Fermi model. Default is 0.25. 

i gap_method: Optimality gap method. Possible values: *none*, ^center^, * corners *, * 

corners . cent er \ Default is ^none\ 

# gap: Accepted optimality gap. Default is 0.05. 

# 

# Returns: 

# Final set of sectors, where the first row yields the solution. 
terminate <— FALSE ; iteration <— 0 

data <— .create. data(targets , receivers) 

sectors <— .create. sectors(data , rho_0 , obj, sensor, b, ’initial’) 
while (! terminate ) { 

iteration <— iteration + 1 

sectors <— sectors[ order(~ sectors$ upper.bound, sec tors $ longest.edge) ,] 
gamma <— sectors [1 ,] 

gamma$ lower . bound <— . lower . bound ( data , gamma, rho_0 , obj, sensor, b, gap_method) 
gamma$ optimality <— gamma$lower . bound / gamma$upper . bound 
if (gammaSlongeSt . edge <= longest_edge I 1 — gamma$optimality <= gap) 
terminate <- TRUE 
else { 

if (verbose) { cat ( pasteO ( ’ I te r at i o n ’, iteration , ':\n’)); print (gamma) ; cat(’\n’)} 

if (gammaSxmax — gamma$xmin <= longest_edge) xs <— c(gamma$xmin, gamma$xmax) 
else xs <— seq (gammaSxmin, gamma$xmax, length, out = 3) 

if (gammaSymax — gammaSymin <= longest_edge) xs <— c(gamma$xmin, gamma$xmax) 
else ys <— seq (gammaSymin, gammaSymax, length, out = 3) 

sectors .prime <— . create . sectors (data , rho_0 , obj, sensor, b, xs = xs, ys = ys, 

name = pasteO ( ’ i t e r ’ , iteration)) 
sectors <— rbind ( sectors [ —1 ,] , sectors . prime ) 


cat ( pasteO (’ Finished after ’, iteration , ’ iterations (evaluated ’, 

nrow (sectors) + iteration — 1, ’ sectors) :\n’)) 


45 





43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 


print (gamma) 
return (sectors) 


create . data <— function (targets , receivers) { 
i Creates the data frame used hy DiBS containing target locations, 

# values and distances to receivers 

M Args : 

i targets: Data frame with object, x, y and value, 
i receivers: Data frame with object, x and y. 

M Returns: 

i Created data frame. 

data <— targets[,c (’ Xyvalue ’) ] 

row . names ( data ) <— targetsSobject 

data[pasteO(’dist.’, as.factor(receivers$object))] <— 
sap ply (row . names ( receivers), function(r) 
supply (row . names (targets), function(t) 

sqrt(( targets[t, ’x’]— receivers[r, ’x’ ])^2+(targets[t, ’y’]— receivers[r, ’y’ ])^2))) 
return (data) 


create . sectors <— function (data , rho_0 , obj , sensor , b, name, xs = NUIX, ys = NULL) { 

# Creates sectors from the given data and evaluates 

M each sector^s longest edge and upper bound. Uses xs and ys as 
i slices. If xs and ys are NULL then it creates the initial split. 

# 

# Args : 

M data: Data set created by function . create. data(). 
i rho_0: Range of the day 
M obj: Objective. 

# sensor: Sensor Model. 

^ b: Diffusivity parameter for Fermi model, 

i name: identifier for the set of sectors 

# xs: Vector of vertical slices. Default is NULL. 

i ys: Vector of horizontal slices. Default is NULL. 

» 

i Returns: 

i Data frame with sector coordinates and upper bounds. 
if ( is . null ( xs )) xs <— seq ( min( data$x) , max(data$x), length. out=3) 
if ( is . null ( ys )) ys <— seq ( min( data$y) , max(data$y), length. out=3) 
sectors <— data.frame( 


sector 

= pasteO (name, 1 

: (( length ( xs ) — 1) * 

( length (ys) 

xmin 

= rep ( xs[— length ( xs )] , 

length ( ys ) — 1) , 


xmax 

= rep ( xs [ - 1 ], 

length ( ys ) — 1) , 


ymin 

= rep (ys[“ length ( ys )] , 

each = length(xs) - 

1), 

ymax 

= rep(ys[-l]. 

each = length(xs) - 

D) 


s e c tor s $ longe s t . edge <— pmax( s e c t o r s $xmax — sectors $xmin , sectorsSymax — s ec tor s $ymin) 
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sectors Supper . bound <— sapply (1: nrow( sectors ) , function(s) { 
switch(obj , avg = value <— 0, min = value <— Inf, 

stop(pasteO(’objective obj, not known ’ ))) 

if (distr == ’cookie’) value <— 0 
for(t in 1: nrow( data )) { 

X <— min(max( sectors $xmin [s], data$x[t]), sectors $xmax [ s ]) 
y <— min(max( sectors Symin [s], data$y[t]), sectors Symax [ s ]) 

Pt <— .Pt(data, t, X, y, rho_0 , sensor, b) 
if(sensor == ’cookie’) { 

value <— value + Pt*data[t, ’value’] 

} else { 

switch (obj , 

avg = value <— value + data[t, ’value’] * Pt / nrow(data), 
min = value <— min(value , data[t, ’value’] * Pt) 

) 

} 

} 

return (value ) 

1) 

return (sectors) 


. Pt <— function (data , t, x, y, rho_0 , sensor, b) { 

# Calculates detection probability Pt for target t with source at position (x, y) 

M Args : 

i data: Data set created by function . create. data(). 
i t: target 

i X: source * s x coordinate 

# y: source ^s y coordinate 
i rho_0: Range of the day 

# sensor: Sensor Model. 

^ b: Diffusivity parameter for Fermi model. 

# 

# Returns: 

# Detection Probability Pt 

d.ts <— sqrt((x — data$x[t])^2 + (y — data$y[t])^2) 

Pt <- 1 

for(r in 4: ncol (data)) { 

rho.tsr <— sqrt(data[t,r] * d.ts) 
switch ( sensor , 

cookie = if (rho.tsr <= rho.O) {Pt <— 0; break}, 

fermi = Pt <— Pt * (1 — 1 / (1 + 10^((( rho . tsr / rho.O) — 1) / b))), 
exp = Pt <- Pt (1 - 10^(-0.30103 * rho.tsr / rho.O)), 

stop(pasteO(’distribution distr, ’" not known ’ )) 

) 

) 

return (1 —Pt) 

) 
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lower. bound <— function (data , gamma, rho_0 , obj , sensor, b, gap_method) { 

# calculates a lower hound for sector gamma 

» 

i Args : 

# data: Data set created hy function . create. data(). 

# gamma: Sector data created in function DiBSO 
^ rho_0: Range of the day. 

M longest_ edge: Termination condition, 

i ohj: Objective. Possible values: *avg\ *min\ 

M sensor: Sensor model. Possible values: ^cookie\ ^fermi\ ^ exp \ 
i b: Diffusivity parameter for Fermi model. 

# gap_method: Optimality gap method. Possible values: ’none*, ’center’, ’corners’, ’ 
corners. center ’. 

# Returns: 

^ Lower bound for sector gamma. 
switch (gap_method , 

none = return (0) , 

center = { 

potentials <— data.frame(x=mean( c(gamma$xmin , gamma$xmax)) , 

y=mean( c (gamma$ymin , gamma$ymax))) 


corners = { 

potentials <— data.frame (x=c (gammaSxmin , gammaSxmax , gammaSxmin , gammaSxmax) , 

y=c (gammaSymin , gammaSymin , gammaSymax , gammaSymax)) 


center.corners = { 

potentials <— data, frame (x=c (gammaSxmin , gammaSxmax, gammaSxmin, gammaSxmax, 

mean( c(gammaSxmin , gammaSxmax))) , 
y=c (gammaSymin, gammaSymin, gammaSymax, gammaSymax, 
mean( c (gammaSymin , gammaSymax)))) 

}, 

stop ( pasteO (’ gap method gap_method , not known’)) 

) 

potentials Svalue <— sapply (1: nrow ( p o t e n t i al s ) , function(p) { 
switch (obj , avg = value <— 0, min = value <— Inf, 

stop(pasteO(’objective obj, not known ’ ))) 

if (distr == ’cookie’) value <— 0 
for(t in 1: nrow( data )) { 

Pt <— .Pt(data, t, p o t e n t i al s $x [ p ] , p o t e nti al s $y [ p ] , rho_0 , sensor, b) 
if(sensor == ’cookie’) { 

value <— value + Pt*data[t, ’value’] 

} else { 

switch (obj , 

avg = value <— value + data[t, ’value’] * Pt / nrow(data), 
min = value <— min(value , data[t, ’value'] * Pt) 

) 

} 
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return (value) 


1 ) 

potentials <— p o te nti al s [ order(— p o te n t i a 1 s $ value ) , ] 
return ( pote nti al s $ value [ 1 ]) 


# Ex amp I e 

# The next examples demonstrates the use of DiBS. 
targets <— data.franie( 


obj ect = c( 

’T1 ’ , 

’T2’ 

, ’T3’ , 

’T4’ , 

’T5 ’ , 

’T6' , 

’T7’ , 

’T8’ , 

’T9’ , 

’TIO’ ) 

X = c ( - 

15.00, 

6.25, 

- 9.25, 

0.75,- 

4.00, 

9.00,- 

9.50 , 

2.75 , 

3.50,- 

4.00) 

y = c(- 

4.75,- 

9.50, 

-11.00, 

8.75 , 

5.50, 

-11.25,- 

15.00, 

11.50, 

-12.25 , 

7.50) 


value = 1) 

receivers <— data.frame( 

obj ect = c ( ’R1 ’ , ’R2 ’ , ’R3 ’ ) , 

X = c( 4.25 , 14.25,- 6.75) , 

y = c(-14.25, 6.25 , 12.50)) 

# 1. fermi model, maximize average detection probability, longest edge <= 0.25 
sectors <— DiBS(targets , receivers , rho_0 = 3, longest_edge = 0.25) 

M 2. cookie cutter model, longest edge <= 0.1 

sectors <— DiBS(targets , receivers , rho_0 = 3, longest_edge = 0.1, sensor= ’cookie’ ) 

# 3. like 1. but longest edge <= 0.001, use center for lower bound 

sectors <— sectors <— DiBS(targets , receivers , rho_0 = 3, longest_edge = 0.001, 

gap_method = ’center’, gap = 0.01) 


DiBS Implementation - This code shows a basic DiBS implementation in R. It features all 
sensor models mentioned in Section 1.2 as well as longest edge and optimality gap termination. 
It assumes that plane rotation has been applied prior to running DiBS. Uneven rectangles are 
handled by stopping sector splitting in the direction that reaches the longest edge termination 
condition. 
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Supplementals 


• Computer Code of the analysis and decision tool - to be used with Appendix A. 

• Computer Code of the DiBS algorithm implemented in R - to be used with Appendix 
B. 

The supplementals are available at Dudley Knox Library of the Naval Postgraduate School 
in Monterey, CA. 
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